INT-PUB 05-14 



in 
o 
o 

p< 



> 

00 



S3 



The Color Glass Condensate and hadron production 

in the forward region 



Adrian Dumitru", Arata Hayashigaki" and Jamal Jalilian-Marian^ 



' a Institut fur Theoretische Physik, J. W. Goethe Universitat, 



Max-von-Laue Strasse 1, 
D- 60 43 8 Frankfurt am Main, Germany 
b Institute for Nuclear Theory, University of Washington, Seattle, WA 98195 

O 

in 
o 

a,: 

Abstract 

" 1—1 , 

^ ■ We consider one loop corrections to single inclusive particle production in parton- 



nucleus scattering at high energies, treating the target nucleus as a Color Glass 
Condensate. We prove by explicit computation that in the leading logQ 2 approxi- 
mation, these corrections lead to collinear factorization and DGLAP evolution of the 
projectile parton distribution and hadron fragmentation functions. In single-inclusive 
cross sections, only two-point functions of Wilson lines in the adjoint and fundamen- 
tal representations (Mueller's dipoles) arise, which can be obtained from the solution 
of the JIMWLK equations. The application of our results to forward-rapidity pro- 
duction shows that, in general, recoil effects are large. Hence, the forward rapidity 
region at RHIC is rather different from the central region at LHC, despite compara- 
ble gluon densities in the target. We show that both the quantum x-evolution of the 
high-density target as well as the DGLAP Q 2 -evolution of the parton distribution 
and fragmentation functions are clearly seen in the BRAHMS data. This provides 
additional strong evidence for the Color Glass Condensate at RHIC. 



1 Introduction 



The rapidity dependence of the recent RHIC data on hadron production in deuteron-gold 
collisions [1] may hint at the emergence of the Color Glass Condensate (CGC) [2] as the 
dominant physics in the forward rapidity region [3]. While the ratio of the hadron trans- 
verse momentum distributions in d + Au versus pp collisions shows a Cronin enhancement 
at midrapidity and moderate transverse momentum [4], the enhancement turns into sup- 
pression as one goes to larger rapidities. The disappearance of the Cronin peak is most 
commonly taken to be due to the quantum evolution (with log 1/x) of gluons in the tar- 
get [5, 6, 7]. Another argument in favor of the CGC formalism is provided by the changing 
centrality dependence of the data as one goes to the forward region. Note also that leading- 
twist next-to-leading order perturbative QCD calculations [8] provide a good description 
of the inclusive distribution at forward rapidities in pp collisions but fail to describe the 
d + Au BRAHMS h~ data by a conventional modification of the leading-twist parton den- 
sities in nuclei (shadowing). The CGC dynamics can be tested further by p + A collisions 
at the LHC, where measurements in both the central and forward rapidity regions should 
be performed. Interestingly, particle production in the forward region of hadron-nucleus 
collisions also plays a key role for the properties of giant air showers from ultra-high energy 
cosmic ray interactions in the atmosphere [9]. Improving our understanding of particle pro- 
duction at large Feynman-x may therefore reveal the type and origin of those super-high 
energy particles. 

Even though the qualitative predictions of the CGC were confirmed by the RHIC data, 
in order to firmly establish the CGC as the cause for the observed suppression of the hadron 
spectra in deuteron-gold collisions and to clarify the role of other scenarios [10] without 
log 1/x resummation, one needs to consider effects in the CGC framework which have been 
neglected so far and which may be significant. One such effect is the recoil of the source 
which radiates the produced gluons. Current approaches to gluon (hadron) production in 
proton-nucleus collisions [11, 12, 13] treat both the projectile proton and the target nucleus 
as a Color Glass Condensate, with the difference that the proton is assumed to be in the 
dilute regime while the target nucleus could be in either the dilute or dense regimes. This 
approach neglects recoil effects since only diagrams that survive in the £ — > limit, with 
£ the momentum fraction of the produced gluon, are included. Gluon radiation and recoil 
effects should also be significant for the production of leading hadrons in deep-inelastic 
scattering from nuclei, where large distortions of the spectrum relative to leading-twist 
calculations have been predicted [14, 15]. 

However, in forward rapidity production of particles, the momentum fraction x of the 
projectile parton is large and treating the proton as a dilute CGC (i.e. subject to BFKL 
evolution [16]) can not be justified. Theoretically, BFKL resummation applies if the pa- 
rameter a s \ogxo/x ~ 1, where xo ~ 0.1 is where the valence degrees of freedom reside. 
In order to gain appreciable longitudinal phase space so that BFKL evolution becomes 
significant, one needs to go much below xq, which is not the case in the forward rapidity 
region, where (by definition) the typical x in the proton wave function is of order xq. It is 
known from HERA data that the saturation scale of a proton reaches ~ 1 GeV only around 
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x ~ 1CT 4 . Since the saturation scale of the proton is so small (of order of Aq C d) at x ~ x , 
even the extended scaling regime is rather small, as well. The production of hadrons with 
p t ^> Aqcd in the forward rapidity region never receives any large contributions from either 
the saturation or the extended scaling regimes of the proton. 

We therefore treat the projectile proton as a collection of quarks and gluons according to 
the parton model and consider scattering of these partons from the target nucleus, which 
is treated as a Color Glass Condensate. In Section 2 we consider one loop corrections 
to quark-nucleus scattering due to gluon radiation and show that recoil effects can be 
large. We then generalize this to include all processes at this order and show that in 
the leading logQ 2 approximation, these lead to DGLAP evolution [17] of the quark and 
gluon distribution functions of the projectile proton, and of their fragmentation functions 
into hadrons. Thus, we prove explicitly that one can treat a high energy proton-nucleus 
collision as scattering of collinearly factorized partons (which evolve according to DGLAP 
evolution equation) in the proton on a dense nucleus treated as a Color Glass Condensate. 

Furthermore, we show that the standard genuinely non-Abelian diagram where the pro- 
duced gluon scatters from the saturated target field gives a contribution which decreases 
by a factor of two toward large rapidity. We then show that QED-like bremsstrahlung, 
which is usually disregarded, contributes about equally at large rapidity and "small" trans- 
verse momentum. Both contributions involve only two-point functions of Wilson lines 
(dipoles [18, 15]), in the adjoint and fundamental representations, respectively. This makes 
it possible to use the RHIC data on hadron production in d + Au collisions to make pre- 
dictions for electromagnetic processes such as photon and dilepton production [4, 19] since 
these electromagnetic processes also involve dipoles in the fundamental representation. The 
evolution of the target wave function with x can then be included by using the solution 
of the JIMWLK equations [2] for the dipole evolution (for first attempts see e.g. [20]) or 
alternatively, by using phenomenological parameterizations of the dipole profile. 

Finally, we apply our results to forward rapidity hadron production in deuteron-gold 
collisions at RHIC energy. Lacking a solution of the JIMWLK equations we resort to 
the phenomenological parametrization of Kharzeev, Kovchegov and Tuchin for the dipole 
profile [6]. We find that the BRAHMS data can be reproduced with a p t - independent 
if-factor. We show that both the DGLAP (Revolution of the distribution and fragmen- 
tation functions as well as the quantum ^-evolution of the high-density target (anomalous 
dimension of its gluon distribution function) are clearly seen in the data. In our opin- 
ion, this strengthens the evidence for the Color Glass Condensate in high-energy d + Au 
collisions [21] substantially. 

Despite the rather good description of the data provided by the above-mentioned for- 
malism, some words of caution may be in order. It has been shown within the NLO pQCD 
framework that at RHIC energy, incoherent leading-twist interactions with partons from 
the nucleus about xa ~ 0.01 contribute significantly [8] (in fact, even overshooting the data 
at forward rapidity). While the contribution from this region can be neglected when very 
small x in the nucleus become accessible at very high energies, at RHIC energy it might 
be necessary to think about additional mechanisms which suppress such interactions [22]. 
In this paper we do not construct any such model, however, but assume the validity of the 
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high-energy approximation inherent in the non-linear x-evolution, namely that the target 
fields at rapidities below that of the probe can be integrated out to yield the effective sat- 
uration scale Q s (y). We return to this point in appendix B, where we show that the 2 — > 1 
like kinematics which arises in the high-energy limit probes very small momentum fractions 
in the target nucleus, on the order of (xa) ~ 10~ 3 for hadrons produced at rapidity Y ~ 3.2 
(BRAHMS) and (x A ) ~ 10" 4 for Y ~ 4 (STAR). 

2 Including recoil in gluon production 




(3) (4) 



Figure 1: Diagrams contributing to gluon production in quark- nucleus scattering. 

The diagrams contributing to gluon production in quark-nucleus scattering are shown 
in Fig. 1. The amplitudes for production of a massless quark with momentum q and a 
gluon with momentum k are given in [23] and read 

M l = -ig 1 ^u{q)i{l 1 + h)l-u{p)t a \V{q t + k t )-{2nf5 2 {q t + k t )} 

M 2 = l g-^u{q) 1 -{t>-h)iu{p)\V{q t + k t )-{2n) 2 5\q t + k t )]t a 

M 3 = tg^u(q)^u(p)d^(p-q)e,(k)t b [U ba (q t + k t )-5 ba (27r) 2 5 2 (q t + k t )] 

[V(q t + l t ) - (2n) 2 5 2 (q t + l t )]t b [U ba (k t - l t ) - 5 ba (2n) 2 5 2 (k t - l t )} , (1) 

where p(q) is the momentum of the incoming (outgoing) quark and z = q~ jp~ the fractional 
light-cone momentum carried by the quark in the final state (£ = 1 — z = k~ /p~ is that 
of the gluon). All other possible diagrams are suppressed by powers of energy and will be 
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ignored. 

V(x t ) = P expfo J dz-A+(x t ,z-)t a ^ , (2) 

U(x t ) = P expfo f dz-A+(x t ,z-)T a ^) (3) 

are Wilson lines, in the fundamental and adjoint representations, respectively, running 
along the light cone and summing up the non-Abelian phases of the colored particles 
propagating through the color-field of the nucleus. To calculate the gluon production cross 
section, we square the sum of the amplitudes (1) and then integrate over the momentum 
of the outgoing quark, q t . 

After squaring the amplitude in eq. (1), it can be shown [24] that some of the terms 
exhibit a collinear singularity upon integration over the transverse momentum q t of the 
final state quark, which arises when the quark and gluon in the final state are collinear 1 . 
Specifically, the terms \Mi\ 2 and |M 3 | 2 are singular while the rest are finite. These two, 
unlike for example |M 4 | 2 , involve only two-point functions of Wilson lines. To leading 
logarithmic accuracy one can ignore the diagrams giving finite terms and keep only the 
contributions of diagrams which are singular, given by 

z(l + z 2 



M 3 | 2 = 16 (p-f 41 \ Z 2) / riV, e«« +h >-" H A (r,) (4) 

Qt J 



with 



H F (r t ) = C F J d 2 bTr c ([V^b - r t /2) - 1] [V(b + r t /2) - 1]) , 

H A (r t ) = lJd 2 bTr c {[U\b-r t /2)-l][U(b + r t /2)-l]) , (5) 

where Tr c denotes the trace over color matrices. 

Integrating over the quark transverse momentum q t , including the phase space factors 
and averaging (summing) over initial (final) state degrees of freedom leads to 



J qA->gX i Q2 . 



H F {ir t ) + H A {r t ) 



(6) 



where Q 2 denotes the factorization scale. In (6) we have introduced the Leading Order 
(LO) quark-gluon splitting function [25] 

£P g/? (0 EE C F [l + (1 - 2 ] • (7) 

It is easy to understand the origin of the two terms in (6); the first term from Fig. 1- 
1 corresponds to QED-like bremsstrahlung, where a free collinear gluon is emitted after 

1 Care should be taken not to confuse collinear and soft singularities. 
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the quark scatters (multiply) from the target. This contribution vanishes in the recoilless 
approximation, i.e. for z — 1 — £ — > 1, and has been neglected in previous computations of 
gluon production in the color glass condensate approach [11, 12, 13]. As we shall show, it 
is important when z ~ £, i.e. for production of gluons with large rapidity or Feynman-x. 

The second term, from Fig. 1-3, corresponds to the case when the incoming quark 
radiates a collinear gluon which in turn scatters (multiply) from the target. This is the only 
surviving contribution for soft radiation, £ — * 0, inherent in the recoilless approximation. 
The effect of recoil for this term is to simply replace the £ — » limit of the LO splitting 
function P g / q by its exact form given in (7), as expected. This leads to a suppression of 
high-energy (forward) radiation by a factor of 2. 

To assess the QED-like contribution ~ H F from a different perspective, note that the 
invariant mass of the intermediate quark propagator in that diagram (see Fig. 1) is given 



This diverges for gluons in the central rapidity region, £ — > 0, and with fixed k t , but 
this diagram vanishes anyways in that limit. On the other hand, for £ ~ z we have 
m 2 ~ (k t — qt) 2 , which is not proportional to energy and becomes small when k t is nearly 
collinear to q t ; this kinematic configuration then gives a contribution that is not suppressed 
by powers of energy. The case £ ^> z, finally, corresponds to "inverse Compton scattering" 
kinematics: a high-energy quark scatters from one or more small-x gluons and emits a high- 
energy gluon which takes over most of the light-cone momentum of the incident quark [26] . 
In order that m 2 not be large (proportional to energy), the quark has to remain nearly 



collinear to the beam: q t <, yz/^q . For collinear gluons with k t ~ £qt/z ~ y £/ ' z q then, 
m 2 ~ 0. 

The collinear logarithm in the first term of (6) can be understood as part of the one 
loop correction to the fragmentation function of the scattered quark while the collinear 
logarithm in the second term corresponds to part of the one loop correction to the gluon 
distribution function. To proceed further, we note that there are pieces in (6) which 
contribute only when the produced gluon has zero transverse momentum. Since we are 
interested in gluons with finite transverse momentum, those pieces can be discarded by 
defining fundamental and adjoint dipole cross sections at impact parameter b as follows 



by 



m2 = tf-z + q 2 t 2q t • k t . 



(8) 





N F (r t ,b) ee j-Tr c (V\b-r t /2)V(b + r t /2)-l) 
N A (r t ,b) = - ] J-_Tr c (U\b-r t /2)U(b + r t /2)-l) 



(9) 



in terms of which the cross section (6) can be rewritten as 



da qA^ 9 x x 



91 [I 

A 2 U 2 



N F (k t /C,b) + N A (k t ,b) 



(10) 



did 2 k t d 2 b (2tt) 
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So far, we have calculated the gluon production cross section from scattering of quarks 
on the target nucleus. Before we can relate this to hadron production in proton-nucleus 
collisions, we need to include all other processes, to the same order in a s , which contribute 
to hadron production. These are, for example, elastic and inelastic production of quarks 
as well as production of gluons from scattering of gluons on the target nucleus 2 . Details 
will be reported elsewhere [24]. Nevertheless, the final result is quite simple and can be 
written down in analogy to eq. (10): 



'dtcPkt d 2 b 
_ da« A ^ x 
'd£d 2 k t d 2 b 

d£d?k t d 2 b 

fagA-.gX 

d£d 2 k t d 2 b 



(2tt) 2 Z P °/*(® 


a s 
2tt 


log 


Q 2 

A 2 


N A (h,b) + 


±N F (k t /Z,b) 


(11) 


( 2vr )2^ P 9/?(0 


a s 
27 


log 


Q 2 

A 2 


N F (k t ,b) + 


±N F (k t /Z,b) 


(12) 


( 2vr )2^ P 9/ 9 (0 


a s 
2tt 


log 


Q 2 

A 2 


N F (k t ,b) + 


±N A (kt/t,b) 


(13) 


(2tt) 2 ^ 9/9 ^ 


a s 
2^ 


log 


Q 2 

A 2 


N A (h,b) + 


^N A (k t /£,b) 


(14) 



with the LO splitting functions as given in [25]. Contributions of processes involving anti- 
quarks are identical to those from quarks to this order. Here, we denote the momentum 
fraction of the produced daughter parton by £, and its transverse momentum by k t . For 
completeness, we also need to include the elastic scattering contributions to quark and 
gluon production [3] given by 



d a iA^ q x 
d^d 2 k t d 2 b 

da gA^gX 

'did 2 kt d 2 b 



(2n) 2 
2 



C6(l-0N F (k t ,b) 
Z6(l-0N A (kt,b) 



(15) 
(16) 



We now reorganize the different pieces from eqs. (11 - 16) in a way which makes the 
connection to the DGLAP evolution of parton distribution and fragmentation functions 
obvious. We consider the diagrams corresponding to the DGLAP evolution of quark dis- 
tribution functions first. They are given by (half of) eq. (15) and the first terms of (12) 
and (13). Multiplying by the bare (parton model) quark and gluon distribution functions 
?o(z/0; 9o(x/Q we obtain 



/ 



I\ q0i l 

f q / P (x/Z,Q 2 )®N F (Z,k t ,b) 



^{ qo ^)[5(l-0 + ^ log ^ P q/q (0 J + go(j)^ log P q/g (0 j N F (£, k t , b) 

(17) 



Q 2 



Q 2 



where f q / p (x,Q 2 ) is the renormalized (DGLAP evolved) quark distribution function in a 
proton. We have used the DGLAP evolution equation for the quark distribution function. 

2 The diagrams considered here correspond to real corrections, one needs to include the virtual correc- 
tions to get the full splitting function. 
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Next, consider evolution of the gluon distribution function. The relevant terms are 
given by (half of) eq. (16), and the first terms of (11) and (14). Putting them together we 
obtain 
i 



,2\ 



Mi-0 + giogf^ /9 (0 



+ %(^log^P g/q (0\N A (t,k t ,b) 



^f 9 / P (x/Z,Q 2 )®N A (Z,k t ,b) (18) 

where f g / p (x, Q 2 ) is the renormalized (DGLAP evolved) gluon distribution function in a 
proton. 

The remaining diagrams lead to DGLAP evolution of the fragmentation functions. We 
start with the quark fragmentation function given by (half of) eq. (15) and the second 
terms in eqs. (11-12) 

k(f ) (V - + | log % P q/q (0) + log % P 9/q (0 



d£ 



A 2 



A 2 



N F (£,k t ,b) 



-> D q (x/Z,Q 2 )®N F (Z,kt,b) (19) 

where D q (x,Q 2 ) is the DGLAP evolved quark-hadron fragmentation function and 
Np{£, k t ,b) = Np(£, k t /C,,b)/^ 2 . The rest of the diagrams give 

/7£ T _ nr. ( m . Cft \ _ t ■ rv fO 2 

N A (Z,k t ,b) 



If 



^(7) (V - + £ log ^ P 9/p (0 ) + ^°(7)S ^g ^ P q/9 (0 



2tt ^° A 2 aiav *') qy C2ir A 2 

-^D g (x/Z,Q 2 )®N A (t,k t ,b) (20) 

where N A (£,k t ,b) = N A (£,k t / 6) / '£ 2 . The final result for one parton radiation can be 
schematically summarized as the following 

-> f q {Q 2 ) ®N F ®D q + q ®N F ® D q {Q 2 ) + / 9 (Q 2 ) ® N A ® D 9 + g ® N A ® D g (Q 2 ) (21) 

where f qtg (Q 2 ), D qt9 (Q 2 ) correspond to DGLAP evolved quark and gluon distribution and 
fragmentation functions, respectively; / , -Do are the bare distribution and fragmentation 
functions, and ® denotes a convolution in x (not shown explicitly). 

The expression (21) is, however, not the complete result as should be clear from the 
presence of bare (non-DGLAP evolved) distribution and fragmentation functions. To get 
the complete one loop result, one needs to consider one additional parton radiation. For 
example, for the diagram shown in Fig. 1-1, one needs to allow for gluon radiation from 
the initial quark line also, which contributes to the one-loop DGLAP evolution of the 
bare quark distribution function. This contribution is calculated in appendix A. The full 
calculation will be shown elsewhere [24], here we just quote the final result for hadron 
production in proton-nucleus collisions 

d(T pA^hX 



dYd 2 Rd 2 b 



(27T) 



7T) 2 J x , 



dx — 

xp Xp 



X 



,X F 



f q/p (x,Q 2 ) N F [—P t ,b] D h / q ( — , Q 2 ) + 

Ob Ob 



f g/p (x,Q 2 ) N A [—P t ,b] D g/h (^,Q 2 ) 

Ob p 1 Ob 



(22) 
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where Y and P t are the rapidity and transverse momentum of the produced hadron while xp 
denotes its Feynman-x. Eq. (22) is our main result. It should be noted that all additional 
parton radiations are taken into account by using the solution of the DGLAP evolution 
equations for the quark and gluon distribution and fragmentation functions. 



3 Forward hadron production in dA collisions at 



We now apply our results to deuteron-gold collisions at RHIC. We use the Leading Order 
CTEQ5 distribution functions for a proton [27] with Q 2 = P 2 , employing isospin symmetry 
to obtain those of the neutron, and the LO KKP fragmentation functions [28] into charged 
hadrons (divided by two). For a deuteron projectile, the difference to fragmentation into 
negatively charged hadrons is small and nearly independent of transverse momentum. 

For our numerical results shown below, we employ the following parameterization 3 for 
the dipole cross section from ref. [6]: 



is the saturation momentum of the nucleus at the rapidity y of the produced parton, 
which can be obtained from its logitudinal momentum x and its transverse momentum 
k t = xP t jxp. The reference point yo = 0.6 specifies where effects due to quantum evolution 
start to become important. Qo is the initial condition for the saturation momentum at yo, 
and the growth rate A ~ 0.3. 

Since data for central collisions is not publicly available yet, we focus on minimum-bias 
collisions. We take an average saturation momentum near midrapidity of Q = 1 GeV [6]. 
The scattering amplitude for a dipole in the fundamental representation, Np, can also be 
parameterized as in (23), with the replacement Q 2 — > Q 2 s Cp/Ca = ^Ql- 

In (23), 7(2/, k t ) denotes the anomalous dimension with saturation boundary condition. 
Ref. [6] employed the following model for the anomalous dimension which reduces to the 
DGLAP anomalous dimension (=1) in the high transverse momentum limit: 



RHIC 



N A (r t ,y) = exp [-I^Q^)]^*)] _ 1 . 



(23) 



Here. 



Qs(y) = Qoexp[\(y - y )/2] , 



(24) 




(25) 



where 



^ h) -(X/2)( y - yo )- 




(26) 



3 Our convention for the sign is opposite to that of ref. [6], which is due to the definitions (9). 
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Figure 2: P t dependence of our results compared to BRAHMS minimum bias data. 



This function vanishes for fixed k t and y — > oo, and so 7 — > 1/2. On the other hand, 
near the boundary to the quantum evolution regime (i.e. at small rapidity y ~ yo such 
that Q s — Qo) the anomalous dimension 7 — > 1, as appropriate for the classical McLerran- 
Venugopalan (MV) model [29]. For both 7 = 1/2 and 7 = 1 the Fourier transform of 
^A{ r t) can De performed analytically: 



N 



N A =1 ( k t) 
7=1/2/7 „2 \ 



A/jr 

d 2 r t e^NT l (r t ) = — exp(-k]/Ql) , 



d 2 r t e i 



ikt-rt j /Y 7=1 / 2 ( 



32tt 



n) 



Ql (1 + 16^)3/2 • 



(27) 
(28) 



Hence, there is an exponential drop with for the MV model [30] with 7 = 1, and a 
power-law decrease if the anomalous dimension 7 = 1/2. Therefore, we expect a steeper 
P t -distribution of hadrons for the classical MV model. In the general case, when 7 depends 
on kt and y, the Fourier transform has to be done numerically. 

In Fig. 2, we show the transverse momentum dependence of the h~ data by 
BRAHMS [31] relative to our results, on a linear scale. The transverse momentum de- 
pendence of the data is reproduced very well. For the correct absolute normalization, we 
need to multiply by a i^-factor of about 1.6; this is not surprising since our calculation 
has been performed at Leading Order (in a s ). In fact, from our point of view it is rather 
comforting that there is room for the expected large NLO corrections. 
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Figure 3: Ratio of the transverse momentum distributions of charged hadrons for Q 2 = 



To test the scale dependence of our LO result, we plot the ratio of the distributions 
for Q 2 /P 2 = 1 and Q 2 /P 2 = 1/4 in Fig. 3. For LO computations there is, in general, a 
monotonic dependence on the hard scale and so "optimal scales" [32], where the result is 
stable against small variations of Q 2 /P 2 , can not be defined. Nevertheless, we observe that 
for Q 2 /P 2 = 1/4 [32] the shape of the distribution does turn out to be nearly the same 
as for Q 2 /P 2 = 1, and that the overall i^-factor drops to ~ 1. This indicates that NLO 
corrections should not change the shape of the distribution by much. In what follows, we 
return to the generic scale Q 2 = P 2 and fix K = 1.6. 

Next, we consider the relative contributions of quarks and gluons in Fig. 4. At rapidity 
~ 3 and p t > 1 GeV quarks clearly dominate [3]. This, of course, is an essential difference 
to the central rapidity region, where gluons contribute more. At yet larger rapidity quarks 
would dominate even at lower p t , and their independent fragmentation should lead to a 
downward shift of baryon- number, which is initially concentrated about beam rapidity [33] . 

In Fig. 5, we compare the full calculation with DGLAP evolution of the distribution 
and fragmentation functions to one where Q 2 has been fixed to 1 GeV 2 . It is obvious 
that DGLAP evolution is important and that it improves the agreement with the data 
significantly. Physically, this is because collinear parton radiation shifts the hard partons 
to smaller momentum fractions, which is the above-mentioned recoil effect, and softens the 
^-distribution of produced hadrons. The effect is clearly seen in the data. This emphasizes 



m/2) 



2 and Q 2 = P 2 . 
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Figure 4: Relative contribution of quarks and gluons from the projectile deuteron. 



that full splitting functions, rather than their soft recoilless limit, have to be employed for 
hadron production in the forward region. 

We also show the result of a calculation within the classical McLerran-Venugopalan 
model, which assumes that the anomalous dimension 7=1. The resulting distribution of 
hadrons is much too soft. It is clear that the data requires the proper quantum evolution 
of the target density, with an anomalous dimension 7 close to 1/2 (with only a weak 
dependence on transverse momentum). This feature is shared by the KKT [6] and IIM [34] 
dipole models, both of which provide a good description of the BRAHMS data. 

Fig. 6 shows the 7r° distribution in minimum-bias d + Au collisions at rapidity y = 
4, which is currently investigated by STAR. Here, x-evolution leads to an increase of 
the saturation momentum by ~ 13% as compared to the BRAHMS kinematics. The 
projectile partons need to carry yet larger momentum fractions (and fragment into faster 
hadrons), hence Q 2 -evolution due to emissions has an even stronger "softening" effect 
on the final hadronic spectra (roughly one order of magnitude at pt = 3 GeV). The Pt- 
distribution at y = 4 is close to exponential for P t : 1 — > 2.5 GeV. At higher transverse 
momentum the convex behavior of the projectile parton distribution and fragmentation 
functions eventually takes over. However, due to the very steep distribution the relative 
deviation of the classical MV model from the full result including quantum evolution in 
x is smaller than for the BRAHMS kinematics. Forthcoming STAR data will test this 
prediction (for preliminary data at intermediate Pt see [35]). 
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Figure 5: Importance of DGLAP evolution of the proton/hadron distribu- 
tion/fragmentation functions and of the anomalous dimension of the target gluon dis- 
tribution. 

4 Summary 

In summary, we have shown that high-energy proton-nucleus collisions can be described as 
scattering of collinearly factorized partons (which evolve according to the DGLAP equa- 
tion) in the proton on a dense nucleus treated as a Color Glass Condensate. We have 
isolated diagrams with collinear singularities which lead to logarithms of Q 2 and proven 
that they satisfy DGLAP evolution. Only two-point functions of Wilson lines, i.e. dipoles, 
contribute to the single-inclusive cross-section; in the future, these could be obtained from 
the JIMWLK small- x evolution equation. 

To apply our results to data from RHIC, we have presently adopted a phenomenological 
dipole parametrization due to Kharzeev, Kovchegov and Tuchin. The minimum-bias data 
from d + Au collisions at forward rapidity obtained by the BRAHMS collaboration can 
be reproduced very well with a transverse momentum independent i\~-factor of 1.6 (for 
q2 _ p2y Hence, NLO corrections are expected to be large but should not distort the 
shape of the transverse momentum distributions obtained from our LO analysis. 

We have shown that in order to reproduce the shape of the measured transverse mo- 
mentum distribution, one needs to account for both (Revolution of the projectile parton 
distribution functions and of the fragmentation functions, as well as for quantum evolution 
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Figure 6: Same as Fig. 5 but for 7r° production at rapidity y = 4. 
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of the small-x field of the target nucleus. Neglect of "recoil" effects from DGLAP evolution 
leads to a significantly harder p t distribution than the data; a purely classical target field 
with anomalous dimension 7 = 1 produces a much softer distribution. The fact that both 
of these key components, and of their interplay, is seen in the data represents important 
evidence for the Color Glass Condensate from RHIC. 

In the future, it would be useful if central rather than minimum-bias data were avail- 
able. For central collisions the saturation momentum of the nucleus increases by ~ 50% 
and so the saturation regime extends to higher transverse momentum. Also, it will be 
interesting to see whether forthcoming data from STAR for ir° production at yet larger 
rapidity, y ~ 4, can be reproduced equally well; we have provided a prediction based on 
the current setup. A better p r resolution will be important to constrain the dipole profiles 
more tightly. Quantitative theoretical ab-initio computations of those universal CGC func- 
tions will hopefully emerge in the near future. We have provided a systematical framework 
which will make it possible to compare them to experimental data. 
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A Radiation from both initial and final state 




Figure 7: One of the diagrams contributing to DGLAP evolution of bare quark distribution 
and fragmentation functions. 



The diagram considered here is shown in Fig. 7. The cross section for this process is 

q - = _J L_ [ ^ f IMP (29) 

y dq-d 2 q t d 2 b 2p- 2(2tt) 3 J (2n)^2ki J (2nf2k^ 11 1 } 

which can be related to hadron production cross section via 

da pA ^ hx f f da qA ^ qX 

Ph dP h -d*P t d% " J dX « qM J dZh D ° iZh) q ~ dq -d V b ■ (30) 

We have defined x q = p~/P~,Zh = P^/q~ and p~,P~ are the momenta of the incoming 
quark and proton respectively, while q~,P^ are the momenta of the outgoing quark and 
hadron. We also define Z\ = {p~ — fcf )/p~ and z 2 = q~ / (q~ + k 2 ). The matrix element is 
given by 

= g 2 u(q)l u t b S Q F (q + k 2 )r f (q + k 2 ,p- k^ip - kj-ft'uip) , (31) 

where 

T f {q,p) = {2n)5{p- - q~) 7 " J d 2 x e i{qt ~ Pt) - Xt [V{x t ) - 1] . (32) 

We need to contract this with the polarization vectors for the two radiated gluons 
e /i(^ij ^i), e u(k 2 , A 2 ) and then square it. Using the polarization tensor J2x e ^(^? A)e l/ (A;, A) = 
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[-9,u + v ^ v H we get 



IMI 



.9 



4 1 



1 



2 (g + k 2 y (p - A^) 4 
Tr D 



r , V^u + VshtJ,^ , Vuh p + r] p k 2l/ 



77 • fci 



77 • A; 2 



Trv 



/Mrf + ^)7 _ (^ " hWtftt - *i)7 _ G* + £ 2 )7 P 
t h t b [V{q t + A; 2t + fc lt - p t ) - (2n) 2 5 2 (q t + k 2t + k lt - p t )\ 



t a t a [V\q t + k 2t + k lt - Pt ) - (2n) 2 5 2 (q t + k 2t + k lt - p t )\ 



(33) 



where Tr D ,Tr c stand for traces over spinor and color matrices. Evaluating the trace of 
the spinors above is algebraically involved but straightforward. Contracted with the gluon 
polarization tensors, it is given by 

[-9,8 + • • •][-&,, + • • •] Tr D [. ■ ■] = 124p'q~(p ■ h)(q ■ k 2 ) }) + + * a) (34) 

^2(1 - - z 2 ) 

where zi, z 2 are defined above. Putting everything together, the partonic cross section can 
be written as 



_ da qA ^ qX g 4 C 



_ g^Up -p- [ dzjl + z\ r dz 2 l+ z\ 
~ (27r) 4g J Zl 1-zJ z\ l-z 2 [Xq 



dq-d 2 q t d 2 b (2tt) z 



Q 



z\z 2 P- 



-) J d 2 x t d 2 y t 



iqr{xt-y t ) 



d 2 k u e iku < xt ~ yt) f d 2 k 2t e ik2t < Xt ~ yt) 



J (2tt) 2 kj t J (27 



) 2 Iht - 



[V{x t )-l\[V\y t )-l\ 



(35) 



Again, the integrations over the transverse momenta k lt , k 2t exhibit collinear singularities 
which lead to logarithms of Q 2 , so that the cross section can be written as 



da qA ^ qX 



1 q f dz\ 



* dq~d 2 q t d 2 b (2n) 2 P-J z x 

S(x q 



^lngc,i±^ 
2tt A 2 1 - z 1 



dz 2 

z 2 



a s Q 2 1 + z 2 

— m —ttCf 

2tt A 2 l-z 2 



Z\Z 2 P z 2 



(36) 



Using (36) in (30) and combining it with the relevant term in (21) gives the quark contri- 
bution to the hadronic cross section 



da pA-+hX 



dyd 2 P t d 2 b 
D (z h ) 



Xp 



dz\ dz 2 dzh . Xp , 

g [ ] 

Z\ Z 2 Zh Z\Z 2 Zh 



^-^ + ^ln^P q/q ( Zl ) 



5(1 -z 2 ) + ^-\nj^P q/q (z 2 ) 



1 ;N F [^,b] 



(2tt) 2 (z 2 z h ) 2 z 2 z h 



(37) 



Defining the DGLAP evolved quark distribution function fq(-^^,Q 2 ) as 



Z 2 Zh J Z-\ 



dz\ Xp 



Zi Z\Z 2 Zh 



S(l-z 1 ) + ^-\n^P q/q (z 1 ) 



(38) 



16 



then changing variables to z' h = z 2 Zh, and defining the DGLAP evolved fragmentation 



function D q {z 



Q 



as 



D q (z' h ,Q 2 )= Do A 

J Zo Zo 



dz- 

leads to the DGLAP evolved hadron production cross section 



5(1 -z 2 ) + ^lnj^P q/q (z 2 ) 



da pA ~ 



■hX 



dyd 2 P t d 2 b 



f dz> h ^rf q ^,Q*)D q (z' h ,Qi 

Jx F Z h Z h 



^z a 



,Pt 



b] 



(39) 



(40) 



°h *h ( 27r ) *h 

It should be noted that one also needs to include diagrams where one integrates over 
the final state quark momenta rather than the gluon momenta as we have done here. 
This would bring in the quark-gluon splitting function P g / q which, combined with P q / q 
above, would then go into the DGLAP evolution of the quark distribution or fragmentation 
function in (38, 39). A further change of variables jz' h in (40) then gives the first 

part of (22). Contribution of the diagrams involving an incoming gluon are similar and 
give the second part of (22). 



B 2 — > 1 kinematics 

In this appendix we elaborate on the 2 — > 1 like kinematics employed above. In particu- 
lar, we show that energy-momentum conservation implies that much smaller momentum 
fractions are probed in the target than for 2 — > 2 kinematics underlying leading-twist per- 
turbative computations (see e.g. [8]). We work in a frame where both projectile and target 
have large light-cone momenta, for example the center of rapidity frame. The various mo- 
menta as used within this appendix are defined in Fig. 8. The incoming parton carries 



p p i 5 




K 

Figure 8: The 2 — ► 1 kinematics. 

momentum p while that of the projectile nucleons is labeled P; the outgoing parton's mo- 
mentum is I while that of the outgoing hadron is q. The momentum of the nucleons from 
the incoming nucleus is K, and k is the exchange between the parton and the nucleus. 
More explicitly 4 , 

P 11 = (P ~ = \[7/2, P + = 0, P t = 0) 

4 We follow the convention from the main text that P and q have large minus-components, i.e. that 
the rapidities of both the projectile and of the produced hadron are large and negative. To comply with 
standard practice, however, in the final figures we quote the modulus of the rapidity of the detected hadron, 
y h = -Y > 0. 
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and 



p v = (p- = I p-p+ = 1 ft = 0) 

z" = (r = P -z+ = J-,z t ) 

^ = (jfe- « 0,A; + = x A K + ,k t ) 

K» = (K- = 0,K + = y/7/2, K t = 0) (41) 
l = p + k (42) 



from energy-momentum conservation. The only approximation made in (41) is setting 
A;~ which corresponds to the eikonal approximation. Using (42), we have k t = l t , 
xP~ = l~ and k + = l + , so that z = x F /x and, finally, 

xa = ^~, (43) 
xps 

where Feynman-x is defined as x F = q~ / ' \J 's/2. This relation can be rewritten in terms 
of the rapidity of the produced hadron, Y (which in the massless limit equals that of its 
parent parton): 

x A = xe' 2yh . (44) 

Here, = —Y > denotes the modulus of the rapidity of the observed final-state hadron. 
In terms of its momentum, its rapidity is given by q~ = q t e Vh / y/2. 

Equation (43) relates x A to the momentum fraction x carried by the impinging projectile 
parton, and to the transverse momentum and Feynman-a; of the produced hadron. Hence, 
we can insert this form inside the integral from eq. (22), then divide by eq. (22) itself, 
to determine the average x A probed in the target nucleus. It is clear from (43) that 
(x A ) -> qf/s ~ 1(T 4 as x F -> 1 (since 1 > x > x F ). The result for both BRAHMS 
and STAR kinematics is shown in Fig. 9. It turns out that for BRAHMS kinematics for 
example, (x A ) ~ 10~ 3 , which is more than an order of magnitude smaller than for the 
2^2 kinematics employed in leading-twist calculations [8]. 

In Fig. 10 we show the kinematic region in x which contributes to the cross section 
given by (22), at hadron transverse momentum q t = 2 GeV and rapidity = 3.2 and 4, 
respectively. Clearly, the cross section is dominated by very large x (labeled in this figure 
as x p ) and very small values of x A . 
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